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Abstract 

The homogeneity of the cosmic microwave background radiation (CBR) is one of the most 
severe constraint for theories of the structure formation in the universe. We investigated 
the effect of the gravitational scattering (lensing) of galaxies, clusters of galaxies, and super- 
clusters on the anisotropy of the CBR by numerical simulations. Although this effect was 
thought to be unimportant, we found that the gravitational scatterings by superclusters 
can significantly reduce the anisotropy of the CBR. We took into account the exponential 
growth of the distance between two rays due to multiple scatterings. The bending angle 
of each ray grows through the random walk process. On the other hand, difference be- 
tween two rays grows exponentially while it is small. This exponential growth is caused by 
coherent scatterings that two rays suffer, and was neglected in the previous studies. The 
gravitational scattering by superclusters reduces the observed temperature anisotropy of 
the CBR at present time approximately by 40-60 % from that at the recombination time 
for angular scale up to a few degrees if the supercluster were formed at z =2-4. 

Subject headings: cosmic microwave background, gravitational lensing, large-scale struc- 
ture of Universe 

1. Introduction 

Observations have shown that the cosmic microwave background radiation (CBR) is 
extremely isotropic. Observed upper limits of the temperature fluctuation, AT/T, of the 
CBR are - 4.5 x 10" 5 at 4'.5 (Uson and Wilkinson 1984), and ~ 2.1 x 10" 5 at V. 15 
(Readhead et al. 1989). Recent observation with COBE showed that the temperature 
fluctuation is (1.1 ±0.18) x 10 -5 at 10 degrees (Smoot et al. 1992), and the observation at 
South Pole showed that the upper bound is 1.4 x 10 -5 at a degree scale (Gaier et al. 1992). If 
the density fluctuation is actually small as suggested by these observations, the structure 
formation models which are consistent with the "observed" temperature anisotropy are 
rather few. For example, the baryon dominant model is ruled out (Peebles and Silk 1990; 
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Gouda, Sasaki and Suto 1989). In addition, the standard CDM model is only marginally 
acceptable (Gouda and Sugiyama 1992). 

During its travel from the last scattering surface, the CBR is gravitationally scattered 
by astronomical objects, such as galaxies, clusters of galaxies, and superclusters. We 
investigated the effect of the gravitational scattering on the anisotropy of the CBR, taking 
into account the exponential growth of the distance between nearby rays through multiple 
scatterings. We found that the temperature anisotropy of the CBR can be reduced by 
40-60% through scatterings by superclusters. 

The smoothing of the CBR by the gravitational lensing has been investigated by many 
researchers (Blanchard and Schneider 1987; Kashlinsky 1988; Tomita 1988; Cole and Efs- 
tathiou 1989; Sasaki 1989; Linder 1990; Tomita and Watanabe 1989; Watanabe and Tomita 
1991). Sasaki (1989) and Linder (1990) gave the mathematical formula of the gravitational 
lensing on the angular correlation function of the temperature anisotropy. 

Kashlinsky (1988) derived equations that describe the multiple gravitational lensing of 
the CBR, and concluded that the original fluctuation was smoothed out on scales up to 
several arcminutes. Cole and Efstathiou (1989) pointed out that Kashlinsky overestimated 
this effect because he modeled galaxy clusters as point masses. They argued that the effect 
is negligible using the formula of increase of the beam width derived by Gunn (1967). 
Blanchard and Schneider (1987) and Watanabe and Tomita (1991) also obtained similar 
results. Tomita and Watanabe (1989) performed numerical simulations of propagation of 
light, and concluded that the effect of the gravitational lensing by clusters of galaxies was 
small. 

Gunn (1967) derived the formula for change of the beam width under an assumption that 
total effect of increase of the beam width is expressed by a superposition of that of small 
scatterings. However, this assumption is not appropriate. If the distance between two rays 
is small, they are scattered coherently. The average increase of the distance between rays 
by one scattering is proportional to its width, since the difference in the deflection angle 
is caused by the tidal force. As a result, the distance increases exponentially by multiple 
scatterings. This exponential growth continues until the scattering becomes incoherent. 
Therefore, evidently, the increase of the distance between rays cannot be expressed by the 
simple superposition. 

In the field of stellar dynamics, the exponential instability of a small difference of initial 
condition has been well known (Miller 1964; Lecar 1968; Sakagami and Gouda 1990; Suto 
1991; Kandrup and Smith 1991; Kandrup, Smith and Willmes 1992; Quinlan and Tremaine 
1992; Huang, Dubinski and Carlberg 1993). However, a clear theoretical understanding 
of this exponential growth was given only recently by Goodman, Heggie and Hut (1993). 



2 



Therefore, it is not surprising that this effect has been neglected in the study of the effect 
of the gravitational scattering on the CBR. 

In this letter, we investigate the exponential growth of the distance between nearby 
rays by numerical calculations, and evaluate its effect on the anisotropy of the CBR. We 
calculated sets of path of two photons in a uniform distribution of the scattering objects. 
Our calculation showed that the angle between two rays increases exponentially up to an 
angle 9 pT given by 9 pT — \/4ivN~2 , where N is the number of the scattering objects within 
the horizon of the present universe in the case of a flat universe. When the angle becomes 
larger than # pr , it increases in proportion to i 1 / 2 . This is because the scattering becomes 
incoherent between rays. 

The angle grows exponentially only when the half- mass radius of the scattering objects, 
rh, is smaller than the projected mean particle distance: d pT — \/4/kN~ % i?n , where i?n 
is the distance to the horizon. We found that superclusters satisfies this condition of the 
exponential growth, rh < d pi , and that 9 pr for supercluster is a few degrees. 

We conclude that the anisotropy of the CBR is smoothed up to scale of a few degrees. 
We evaluated the angular correlation function using the formula derived by Wilson and Silk 
(1981), and Sasaki (1989). We found that the gravitational scattering due to superclusters 
can decrease the anisotropy of the CBR approximately by 40-60% if the beam width of 
antenna is not much smaller than the intrinsic angular scale of the fluctuation. 

2. Numerical Calculations 

We calculated sets of paths of two photons in a uniform and random distribution of 
the scattering objects. We used the post-Newtonian equations of motion for photon (Will 
1981): 

^ = 2V0-4n(n-V0), (1) 

where x is the position of photon, <p is the gravitational potential, and n is the initial 
direction vector which satisfies |n| = 1. This post-Newtonian approximation is always 
valid in a real universe, unless we consider a universe dominated by supermassive black 
holes. The maximum deflection angle of a scattering by an object is determined by the 
depth of its potential well. The depth of the potential well of galaxies, clusters of galaxies 
or superclusters is by far smaller than c 2 , where c is the light velocity. In the numerical 
calculations, we ignored the second term of the right-hand side of equation (1), since our 
purpose is to obtain the trajectory of a photon. The second term changes the deflection 
angle through the change of the light velocity. Since the changes of the velocity and 
the deflection are both 0((p/c 2 ), the contribution of the second term to the deflection is 
0[((p/c 2 ) 2 ], which is negligible. 
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We calculated the trajectory of a photon using the gravitational potential expressed as 



where r SJ - = x SJ - — x, G is the gravitational constant, N and M are the number and 
the total mass of the scattering objects, x SJ and m s j are the position and mass of the 
j-th scattering object, R is the radius of the sphere in which the scattering objects are 
uniformly distributed, and M T j(r) is the mass contained within a sphere of radius r and 
Af r j(r = oo) = m s j. We used the system of units in which G — R — c — 1 and M — 1/2. 
If we specify this sphere as a flat universe, the unit of time is y/2/Ho and the radius R is 
V^c/iJg, where Hq is the Hubble constant. The masses of the scattering objects were set 
to be equal to each other, i.e., m s j — M/N. The second term of the right-hand side of 
equation (2) was introduced to cancel the global harmonic potential due to the scattering 
objects. We started the numerical integration of the trajectory of photons at the surface 
of the sphere. The initial velocity vector, Vq, of a photon points to the center of the 
sphere and |vo| = 1. To calculate the growth of the distance between initially nearby rays, 
we calculated the trajectory of two photons from slightly different initial conditions. We 
performed calculations for three cases: N = 128, 1024, 8192. We used GRAPE-2A (Ito 
et al. 1993), a special-purpose computer for iV-body problem, for the calculation of the 
gravitational force. We integrated the orbits with the 4-th order Runge-Kutta scheme with 
an automatic timestep adjustment (Press et al. 1986) and the Hermite scheme (Makino 
and Aarseth 1992). The details of our calculation are presented elsewhere (Fukushige et 
al. 1994). 

In figure 1, growth factor, «i s = #1.5 /#o, is plotted against 9q/9 pt where angles #0 an d 
#1.5 are the median values of 200 orbit pairs of 9 at t = and 1.5, respectively. In these 
calculations, the scattering objects were set to be point mass (i.e. M T j(r > 0) = m s,j)- 
Figure 1 indicates that the behavior of cni.5 is divided into three cases, depending on the 
value of 9o/9 pi , i.e., exponential, transitional, and diffusive region. If 9q/9 pi < 0.1, the 
growth factor, a, is about 15-20, independent of the initial condition. In this region, 
the angle grows exponentially. If 0.1 < 9q/9 pt < 1, the growth factor cti.5 is roughly 
represented by 9 pi /9q. In this region, 9 increases up to 9 pi . If 9q/9 pt > 1, the growth 
factor, «i.5, is of the order of unity. These behaviors are independent of the number of the 
scattering objects, N. 

The angle grows exponentially only when the half-mass radius of the scattering objects, 
rh, is smaller than the projected mean particle distance, d pT . We calculated the orbits with 
different mass distributions for scattering objects in order to investigate the effect of the 
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size. In figure 2 the growth factor is plotted against r^/d pi for the King model (King 1966) 
with dimensionless central potential Wq — 5, 8, 12 (where the concentration parameter are 
log[r t /r c ] = 1.0, 1.8, 2.7 for Wq — 5, 8, 12, respectively, where r t is tidal radius and r c 
is core radius) and the Plummer model (Plummer 1911). The number of particle, N, of 
scattering objects was 1024 and the initial angle was chosen so that 9q/9 pt <C 1. The 
behavior is almost independent of the softening law. The slope of decrease of the growth 
factor for larger softening depends on the mass distribution, and is a little smaller for more 
concentrated distribution. 

3. Smoothing of Temperature Anisotropy of the CBR 

In this section, firstly, we derive a formula for the growth factor, a, in the expanding 
universe. Next, we estimate the anisotropy of the CBR in term of the angular correlation 
function. Here, we assume that the universe is flat (fi=l) for simplicity. 

The growth factor, ckls 5 of the angle between two photons that come from the last 
scattering surface is calculated as «ls = [w(zf) + 9(zf)Rfi,]/(9qRqi i ), where w and 9 are 
the distance and angle between photons at given redshift z, respectively, 9q — 9(z — 0), zp 
is the redshift when the scattering object were formed, and i?oL an d Rfl are the distances 
to the last scattering surface from z — and z — z-p, respectively. The relation among the 
above parameters are illustrated in figure 3. 

The distance w and angle 9 are approximately given by the equations: 

dw w . dw 

— = -ce --, w\ z=o =0, and c9 = --. (3) 

The e-folding time, t e , is given by t e — r/^/Gp (Goodman, Heggie and Hut 1993), where 
p is the mass density of the scattering object. Using the result of our calculations, r is 
given by t — 3v / 3(4v / 27rlnai.5) _1 . For the case that the scattering object is point mass, 
r ~ 0.18. Solving the equations (3), we obtain the growth factor «ls : 

a LS = (l + z F )i^, (4) 

where rj — tt£ (67r) - 2T -1 , and Q s is the density parameter of the scattering object. 

In Table 1, we summarize estimates of rh/c? pr and 9 pi for galaxies, clusters of galax- 
ies, and superclusters. The projected mean particle distance, d pi , is calculated as 

i — - — - 

d pT — (3rf sep ) 2 i? H 2 f2 s 4 ^sep- The time scale of the exponential growth, t e , is scaled by 
crossing time, t CT . Since the radius of the sphere to which the scattering objects are 
projected is ~ ct cr , it becomes a function of the density parameter, fi s . Therefore, the 
distance d pT is also a function of fi s . The details are discussed in Goodman, Heggie, and 
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Hut (1993). We adopt i?n(= 2c/3i?o) of 2h 1 Gpc, where h is the Hubble constant in unit 
of 100kms _1 /Mpc. 

Using figure 3, we estimate the growth factor, «ls, and show then in Table 1. The 
gravitational scatterings by superclusters increase the angle by a large factor. Here, we 
assume that the superclusters have a large fraction of the mass in the universe. This 
can be justified from the observation showing that the universe consists of the void with 
very low density and the structures of the scale of supercluster. We adopted the size of 
supercluster of 10/z -1 Mpc. The superclusters typically have elongated shapes. The effect 
of the gravitational scattering is determined by the length of the shortest axis. It is unclear 
whether galaxies (fi s ~ 0.2) and clusters of galaxies (fi s ~0.1) have significant effect or 
not. 

In the following, we give a quantitative estimate of the anisotropy of the CBR in terms 
of the angular correlation function using the formula derived by Wilson and Silk (1981) 
and Sasaki (1989). Temperature fluctuation of the CBR averaged over the beam pattern 
of an antenna is given by 
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(9- a) 



2[C(0-a)-C(9-a)i 



(5) 



where 



C(9;a)~ 



1 



2a s 



ipC{ip) exp 



1 



(e 2 + v 2 ) 



In 



2a s 



dip, 



and a is the beam width of the antenna and Iq is the 0th order modified Bessel function. 
The angular correlation function, C(9), is calculated as 



C(9) = C(0) 



1 + 



G(9)9 2 



21-5 



exp 



9 2 



I 2(02 + G{9)9 2 ) 



(6) 



where G(9) — [«ls(#) — I] 2 - Here, 9 C represents the coherence angle of the intrinsic temper- 
ature fluctuation determined by e.g. the Silk damping (Silk 1968). We use an interpolation 
formula: G{9) — 9^ T /[9 2 + 9 2 T G(0)~ 1 ]. This formula well expresses the behavior of G{9) in 
the entire region of 9, as seen in figure 1. In figure 4, the average temperature anisotropy is 
plotted against 9 for the cases of 9 C — 4', a — 1', zp — 2, 4 and 10. For the case of zp —2-4, 
the observed temperature anisotropy is smaller than the intrinsic one approximately by 
40-60%. 

We thank Tomoyoshi Ito, who developed GRAPE-2A, and Shaun Cole for many critical 
comments on the original manuscript. This research is partially supported by the Grand- 
in-Aid for Specially Promoted Research (04102002) of The Ministry of Education, Science, 
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Table 1. The growth factor of angle between photons coming from the last scattering surface 



object 




^sep 


a 




nJd pT 


#pr 


r 


V 






(h~ 


_1 Mpc) 


(/i^Mpc) 


(h- 


^Mpc) 











z F = 2] 


[z F = 4] 


galaxy 


0.02 


5 


0.2 


0.65 


0.031 


0.50 


0.18 


0.56 


2.5 


3.9 


cluster of gal. 


1 


35 


0.1 


14 


0.071 


7.8 


0.19 


0.38 


1.9 


2.5 


supercluster 


5 


150 


1.0 


71 


0.070 


122 


0.19 


1.20 


7.3 


18.3 



Figure Caption 

Fig.l Growth factor, a.\^ = Oi.s/Oq-, plotted against #o/$pr- The triangles, squares and 
pentagons are for N — 128, 1024, and 8192, respectively. 

Fig. 2 Growth factor, «i.5, plotted against r^/dpj- for different mass distribution of scatter- 
ing objects. Solid curves indicate the King models. The triangles, squares and pentagons 
are for Wq — 5, 8, and 12, respectively. Dotted curve indicates the Plummer model. 

Fig. 3 Growth of angle between photons coming from the last scattering surface. The 
redshift zls is that of the last (Thomson) scattering surface, zy is the redshift when the 
scattering object were formed, w and 6 are the distance and angle between two rays, 
respectively. 

Fig. 4 Averaged temperature fluctuation plotted against the chopping angle, #, for 9 C — 
4' and a — V . The thin curve indicates the intrinsic temperature fluctuation with no 
gravitational scattering. The thick curves indicate the observed temperature fluctuation 
for zy — 2,4 and 10. 
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